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ABSTRACT 

We present the dust properties of a small sample of Virgo cluster dwarf galaxies drawn from the science demonstration phase data set of the 
Herschel Virgo Cluster Survey. These galaxies have low metallicities (7.8 < 12 + log(0/H) < 8.3) and star-formation rates < 10" 1 M yr~'. We 
measure the spectral energy distribution (SED) from 100 to 500 /im and derive dust temperatures and dust masses. The SEDs are fitted by a cool 
component of temperature T < 20 K, implying dust masses around 1 s M Q and dust-to-gas ratios D within the range 10 3 -10 2 . The completion of 
the full survey will yield a larger set of galaxies, which will provide more stringent constraints on the dust content of star-forming dwarf galaxies. 
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Late-type dwarf galaxies are metal-poor systems, and they rep- 
resent unique environments to investigate the properties of dust 
in the low-metallicity regime. Mid/far-infrared (MIR/FIR) and 
submillimetre (submm) observations have provided information 
about the different dust components in dwarfs, showing that 
even metal-poor galaxies may host a significant amount of dus t 
(iThuan et alJll999t iHouck et alJ l2004: Eng elbracht et~al"1l2008l) . 
While the majority of the stud ies to da te have focused mainly 
on bright and isolated dwarfs dGalliano et a l. 2003; Hun t et alj 
l2005tlGaUiano et alJl2005blGalametz et alJ2009h . little is known 
about the interplay between the environment and the dust prop- 
erties of low-mass and low-metallicity systems. 

It is still unclear whether most of the dust mass in these 
systems is at very low temperature (T < 10 K), as the excess 
emission in the submm obs erved in some dwarfs m ight sug- 
gest (Madden 2002; Ga lliano et all2005t iGalametz e t al. 2009). 
This cold dust component may either be associated with the 
star-forming regions , residing in clumpy molecular complexes 
(Gall iano et alJl2003h . or extend beyond the optic al disc follow- 
ing th e distribution of the neutral hydrogen (Hi; iPopescu et alJ 
2002). However, the main objection to this interpretation is the 
high dust mass associated with the cold dust, and alternative ex- 
planations have been propos ed (seelBendo et alj2 010, and refer- 
ences therei n). For example, Lisenfeld et alj d2002l) used a mixed 
dust model dDesert et alJ 19901 of large and small grains, finding 
that a high abundance of small grains (compared to the Milky 
Way) with a shallow emissivity (J3 < 2) may ac count for the ex- 
cess submm emission in NGC 1569. IZhu et alJ (J2009) applied a 
similar approach to the starburst spiral galaxy NGC 3310 with- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



out introducing a cold dust component. Only additional mea- 
surements in the FIR/submm regime and the analysis of a larger 
sample of galaxies can help us properly assess the issue of dust 

at low temperatures. 

The Herschel Space Observatory (Pilbratt et al.l 120101) al- 
lows us for the first time to study the wavelength range be- 
tween 200 and 500 fim, where the emission from cold dust is 
expected to be the most prominent. As part of the Herschel 
Open Tim e Key Project, the Herschel Virgo Cluster Survejjj 
(HeViCS; iDavies et all 12010) will map an area of 64 square 
degrees of t he Virgo cluster w ith PACS (Poglit sch et al.ll2010h 
and SPIRE (iGriffin et al.ll2010h to investigate the dust content of 
the different morphological types within the cluster. The dwarf 
galaxy population of the Virgo cluster is dominated by dwarf 
ellipticals, but also contains a non-negligible fraction (~ 10%) 
of late-type dwarfs w ith signs of current star formation activity 
(Binggeli et al.l ll987l) . In this Letter, we present the dust prop- 
erties of three of these galaxies detected with PACS and SPIRE 
in the HeViCS Science Demonstration Phase (SDP) data set, a 
4° x 4° field covering the central region of the cluster. 



2. Sample selection 

The HeViCS SDP field contains 28 late-type low-luminosity 
galaxies classified in the Virgo Cluster Catalog (VCC, 
Binggeli et al.lll985l) as Sm (7), Im (12), blue compact dwarfs 
(BCDs) (6), and dwarf irregulars (3). Three of the BCDs have 
been detected with PACS and SPIRE and are analysed in this 
work: VCC 562, VCC 1 179, and VCC 1356. They have B mag- 
nitudes between 15.5 and 16.7 (Mb between -14.5 and -15.6), 
and Hi masse s of 10 7 - 10 8 M Q assuming a distance to Virgo 
of 16.5 Mpc (Mei et al. 2007). We inferred star- formation rates 
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Fig.l. SDSS and Herschel images of VCC562, VCC 1179, and VCC 1356 at 100, 160, 250, 350, and 500 urn. The field size is 
150". The FWHM of the beam at each wavelength is displayed at the bottom-right corner. The different shape of the PACS beams 
depend on the single-scan (ellipse) or crossed-link scan (circle) coverage of the area. 



(SFRs) of between 10 _1 an d 10~ 2 M^yr" 1 from extinction- 
corre cted Ha and UV fluxes dGavazzi et alJ 2004; Boselli et alj 
2009). Althou gh these dwarfs are compact, thus classified as 
BCDs by Bin ggeli et all (1 19851). their SFRs are lower than the 
typica l values for BCDs dHunter & Elmegreenl2004UZitrin et alj 
2009). Metallicities were estimated using th e Sloan Digitised 
Sky Survey (SDSS) spectroscopic database (lYork et al. 2000). 
Since the [On] /13727 nebular line is outside the observed wave- 
length range of the SDSS spectra, we estimated the oxyge n 
abundances following the method of Pilyugin & Thuanl (|2007). 
Our estimates agree, within the errors, with the measure- 
ments available in the literature f or VCC 56^3 and VCC 1 179 
dVflchez & Iglesias-Paramd l2003). The basic parameters of the 
three galaxies are given in Table 1. 

Among the Sm and Im types, only one bright IBm galaxy, 
IC3583 (VCC 1686), was detected in the SDP field, but given its 
different properties from the three dwarfs, it will be analysed in 
more detail in a separate work. 



3. Photometry 

The SDP data were re duced using the Level I procedures de- 
scribed in Pohlen et al. (2010) . Mor e details on the observations 
can be found in Da vies etaTI (l2010h . The angular resolution for 
PACS in fast scan parallel mode is 7"xl2'/7 and ll'.'6xl5'.'7, at 
100, and 160 /urn, respectively. For SPIRE, the PSF FWHM is 
18'.'1, 25'.'2, and 3679 at 250, 350, and 500 /mi respectively. 
Pixel sizes are 372, 674, 6", 10", and 14" for the five bands. 
The rms levels of empty sky regions at 100, 160, 250, 350, 
and 500 //m are ~2mJypix _1 , ~5mJypix~ 1 , ~8mJybeam _1 , 
~7mJybeairT 1 , and ~9mJybearrr 1 , respectively. 



2 Note that according to Vilchez & Iglesias ^Paramol J2003) the un- 
certainties in the measurement of VCC 562 are also compatible with a 
higher metal abundance, 12 + log(0/H)=8.4-8.6. 



FigureQ]compares PACS and SPIRE data of the three dwarfs 
to the SDSS images. VCC 1179 is not detected at 100 /mi, nor 
VCC 562 at 500 /mi. The other galaxies are marginally detected 
at 500 /mi with a low signal-to-noise ratio (S/N) (4cr and 5cr for 
VCC 1 179 and VCC 1356, respectively). The VCC 1 179 images 
show an additional feature to the north extending beyond the 
edge of the optical disc. It is not clear whether this feature is 
associated with the galaxy or not, but it also does not seem to be 
related to any background source in the SDSS images. 

Photometry from 100 to 500 /mi was derived by means 
of the standard growth curve analysis. The aperture is cen- 
tred on the brightness peak, and the sky level is determined 
by considering the value that minimizes the radial variation 
in successively larger apertures. At large radii, the photometry 
should be relatively constant, and the asymptotic value corre- 
sponds to the total flux. The value of the sky thus derived is 
always within 1.5 standard deviations of the sky background 
measured in empty regions around the galax ies. The uncertainty 
in the flux density is assumed to be 30% (Boselli et alj|2010l : 
ISwinvard. Ade. Baluteau. et alJl2010h . 



4. Dust temperatures and masses 

To derive the temperatures of the dust, we fitted the 
PACS+SPIRE SEDs with a single modified Planck function and 
an emissivity law k v oc y 2 (Fig. 2, upper panels). We constrained 
the dust masses using the 250 //m flux densities, and the emis- 
sivity k v = 4.67 cm 2 g" 1 at 250 /mi dLi & Draine 2001). The 
resulting dust masses are approximately 10 5 M , and the dust 
temperatures are around 20 K (see Table 1). 

However, Fig. 2 shows that 500 /mi fluxes of VCC 1 179 and 
VCC 1356 tend to be underestimated by the single-temperature 
fits. This difference cou ld be due to e ither thermal or non- 
thermal radio emission dCondonl 1 1992k iHunt et al] l2005h, an 
additional cold (~10K) dust component (Galli ano et al J 12003: 
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Table 1. The properties of the three star-forming dwarf galaxies in the Virgo cluster under study. Dust temperature, masses, and 
dust-to-gas ratios were obtained by fitting a single modified black-body function with/? = 2. The distance to Virgo is assumed to be 
16.5 Mpc. 



ID 


mag 


12+log(0/H) 


SFR 

M G yr' 


log(M„ ; /M ) 


T \BB 

K 


logOVWMo) 


log(£> ) 


VCC 562 
VCC 1179 
VCC 1356 


16.74 
15.46 
16.15 


7.8+0.2 
8.3+0.2 
8.0+0.3 


4.6 xl0~ 2 
9.4 xl0~ 2 
4.6 xl0~ 2 


7.67+0.03 
7.39+0.06 
8.35+0.01 


20.2+0.2 
16.2+0.6 
17.9+1.1 


5.18+0.13 
5.48+0.14 
5.53+0.16 


-2.49+0.13 
-1.91+0.15 
-2.82+0.16 



Galametz et al. 2009), an enhanced abundance of small grains 
dLisenfeld et al.ll2002i) . or the different optical properties of the 
amorphous dust grains dMenv et al.ll2007l) . 

We first considered the possibility of radio emission. 
Following Condon (1992), we estimated the non-thermal and 
thermal radio flux at 500 /mi that would be expected given the 
SFR of each galaxy. Because of the steep frequency fall-off 
of the non-thermal radio component toward higher frequencies 



(- 



-0.8 



), the estimated thermal flux at 500 //m is higher than the 



non-thermal one by a factor of ~9. Nevertheless, in both cases, 
given the low SFR of these BCDs (< 10" 1 M yr _1 ), the expected 
radio flux at 500 /mi is of the order of 200 //Jy, and thus cannot 
be responsible for the 500 //m excess. 

The emission from an additional cold dust component, 
~10K, would peak between 200 and 300 //m and could cause 
the excess we observe at 500 /mi. To explore this possibility, 
we refined the fiducial SED model by adding a second modi- 
fied black body at a lower temperature. For this two-component 
grey -body fit, we also used/? = 2. Since we have four points in 
the SED of VCC 1 179 and the uncertainties in the flux densities 
are large, the fit could not be tightly constrained, and we discuss 
only the results obtained for VCC 1356. More than one set of 
parameters were able to provide a reasonable fit to the SED of 
this galaxy. We inferred mean temperatures of T c = 8 K and T„. 
= 19. 1 K for the two dust components and the fit for these values 
is displayed in Fig. 2 (lower panel). The corresponding total dust 
mass is Mj = 1.6 x 10 7 M . This high mass is inconsistent with 
current chemical evolution models (see next section). The least 
extreme value of cold dust temperatures providing a reasonable 
fit (T c =11 K, with T u , = 19.6 K) still yields a total dust mass 
Mj — 2 x 10 6 M , which is six times higher than the value for 
single-temperature fit. Fitting the data with a more sophisticated 
dust model to test alternative explanations is beyond the scope of 
this Letter; this analysis will be performed in a future work when 
the observations of the field are completed. At this stage, given 
the low S/N of the detections at 500 fim and the large error bars, 
it is difficult to discriminate between the different scenarios. 

5. Dust-to-gas mass ratios 

The dust-to-gas-mass ratio D (Mdust/M gas ) not only provides 
information about the amount of metals that are locked in 
dust grains, but also provides an indication of the star- 
formation history of a galaxy, reflecting the n et balance be- 
tween the formation and destruction of dust (Hirashitaet al. 
2002). The rati o D is known to c orrela t e with the oxy- 



gen abundanc e (Lisenfeld& Ferrara 1998; Edmunds 2001; 
iHirashita et al.l f2002: Ja mes et alJ20 02). Here we check whether 
the gas and dust masses we obtained for these galaxies in the 
Virgo cluster are compatible with the predictions of dust forma- 
tion models. 

To derive £), we used n eutral hydrogen masses available 
from the ALFALFA catalog dGiovanelli et alJl2007t iKent et alJ 



2008, see Table 1). The angular resolution of the 21 -cm data 
(3'.3 x 3'.8) includes all dust emission detected in these objects. 
We assume the total gas mass is given by the atomic component 
only: independent of the lack of CO detections, this assump- 
tion is justified by the low FL-t o-Hi mass fractions expected 
in metal-poor enviro nments (e.g., Roberts on & Kravtsovl [2008: 
iKrumholz et aT]|2009albl) . 

The resulting values of T> are reported in Table 1 . Figure [3] 
shows D for the Virgo BCDs (filled circles) as a function of neb- 
ular oxygen abunda nce. Also illustrated in Fig. [3] are data from 
the li terature (Jame s et al J 120021; iHuntet al J 1200 5; Walte r et al] 
2007). Selected model predictions for t he T) -metallicity cor - 
relati on are also sho wn (Edmunds 2001t IHirashita et aD l2002). 
The Hirashita et al. (2002) models, shown as long-dashed and 
dot-dashed curves, depend on the dust destruction efficiency, 
which begins to be effective when the oxygen abundance 12 + 
log(0/H) ~ 8. Were the interstellar dust mass an approxi- 
mately c onstant fraction o f th e ISM metal abundance, as pro- 
posed by Edmunds (20011) and iJames et al.l d2002h . the relation 
between D and O/H would be linear. A solid line indicates their 
prediction for dust production by SNe, and the short-dashed 
line represents that by only evolved low- and intermediate- 
mass stars. Alternatively, a non-linear trend is predicted by 
iLisenfeld & Ferraral d 19981) and IHirashita et all d2002l) because 
of the additional effects of outflows and dust destruction effi- 
ciency. The two most metal-poor objects plotted, IZwl8 and 
SBS 0335-052 (at 12 + log(0/H ) -1.2), have a v ery low dust- 
to-gas ratio for their metallicity (Huntetal. 2005), but are ap- 
proximately consistent with linearity. 

Our estimates of D for the three dwarfs, based on the single- 
temperature fit of the data, are indeed consistent with what 
would be expected by the linear models. If these galaxies hosted 
a cold dust component with T c < 1 1 K, as discussed in the pre- 
vious section (and shown in Fig. [3]), the corresponding D would 
be much larger than the model predictions. In this case, a signifi- 
cant missing gaseous component (H2, cold Hi), a few times more 
massive than the amount of detected Hi, should be introduced to 
explain this discrepancy. 



6. Comparison with far-IR properties of known 
star-forming dwarfs 

The brightest and most metal-poor BCDs, such as IZwl8 and 
SBS 0335-052E, have been studied with Spitzer, and dust 
has been detected e ven in these l ow abundance environments 
dHouck et al] 12004 IWu et all 120071) . While the SED peak of 
SBS 0335-052E is shifted to shorter wavelengths (A ~ 28 /mi; 
iHouck etaL 20041). I Zw 18 has also been detected at 70 /mi 
dWu et al.ll2007t lEngelbracht et al.ll2008l) . suggesting that a cool 
dust component might be present in this galaxy. Hunt et al. 
(2005) studied the global SEDs of seven BCDs showing that far- 
IR emission in these galaxies peaks at or shortward of ~60 //m. 
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Fig. 2. SEDs of the VCC dwarfs compared to a single modi- 
fied black body with emissivity index /J — 2 (top), and a two- 
component model fit with (3-2, T w = 19.1 K, T c = 8 K for 
VCC 1356 (bottom). 

However, the dwarfs detected with Herschel are more 
"quiescent" than t he typical BCDs analysed in these studies. 
Walt er et al.l (|2007) observed a sample of dwarf irregular galax- 
ies in the M 8 1 group with SFRs and metallicities more similar 
to the Virgo dwarfs. They derived dust masses and temperatures 
of the M 8 1 dwarfs using only the 70 and 160 /jm Spitzer/MlPS 
bands. Their dust-to-gas mass ratios, taking into account the to- 
tal Hi mass of the galaxies, are displayed in Fig. \3\(plus signs), 
and are about one order of magnitude smaller than those of the 
Virgo dwarfs with a similar metal abundance. Final ly, six Im- 
BCD galaxies in Virgo analysed with ISO data by Popescu et al. 
(2002) show similar properties to the dwarfs in the current work. 
In particular, their 60-170 fim SEDs indicate the presence of 
warm dust and a cooler component with a median temperature 
~18 K, in agreement with our results. The resulting dust-to-gas 
mass ratios in a few cases are rather large, ^0. 1 . 

7. Conclusions 

We have presented PACS and SPIRE observations of three star- 
forming dwarf galaxies in Virgo detected in the far-IR/submm 
regime with the science demonstration phase data set for the 
HeViCS survey. The data indicate the presence of cool dust with 
a temperature S20 K and dust masses around 10 5 M Q . We have 
discussed the possibility that these galaxies host an additional 
cold (T < 10 K) component to explain the excess at 500 //m in 
two of the dwarfs. However, the low S/N of the 500 fim detec- 
tions precludes us from drawing firm conclusions. The comple- 
tion of the full area of the survey will enable us to place more 
stringent constraints on the dust content of star-forming dwarf 
galaxies in a dense cluster environment. 
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